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funded  project  titled  "Subsurface  Signature  of  the  Internal  Wave  Field  Radiated  by 
Submerged  High  Reynolds  Number  Stratified  Wakes"  (Award  Number:  N00014-08-1- 
0235).  The  enclosed  document  is  a  revised  version  of  the  original  report  sent  on  May 
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Technical  Objectives 

The  research  supported  by  ONR  grant  number  N00014-08-1-0235  was  conducted  in  two  phases,  as 

indicated  above.  The  primary  objectives  of  each  phase  are  listed  below. 

STAGE-1,  “Internal  Waves  Generated  by  High  Reynolds  Number  Stratified  Wakes”: 

•  Characterize  the  internal  wave  field  emitted  by  a  stratified  turbulent  wake  during  the  Non- 
Equilibrium  (NEQ)  regime  and  the  particular  modifications  associated  with  high  Reynolds  numbers 
and  the  secondary  instabilities  therein. 

•  Determine  the  remote  signature  of  the  wake-radiated  waves,  as  modeled  by  a  2-D  internal  wave 
beam,  by  examining  their  interactions  with  variations  in  the  background  stratification  profile. 

STAGE-2,  “Subsurface  Signature  of  the  Internal  Wave  Field  Radiated  by  Submerged  High  Reynolds 

Number  Stratified  Wakes”: 

•  Characterize  the  low-wind-condition  subsurface  signature  of  a  single  high-frequency  internal  wave 
beam(s)  in  a  uniform  or  vertically  variable  stratification,  the  latter  typical  of  littoral  ocean  conditions. 

•  Correlate  qualitatively  and  quantitatively,  over  a  range  of  Reynolds  and  Froude  numbers,  the  patterns 
of  subsurface  motion  observed  above  a  stratified  turbulent  wake  with  the  flow  patterns  inside  the 
wake  core  to  infer  the  state  (“age”)  of  wake  turbulence. 


Results 

The  results  obtained  in  pursuit  of  the  above  objectives  may  be  grouped  in  the  following  categories,  on 
each  of  which  we  further  elaborate  below: 


1 


lb 


(i)  Mean  flow  generation  and  mitigation  in  large-amplitude  internal  wave  packets. 

(ii)  Interaction  of  an  internal  wave  beam  with  a  model  oceanic  pycnocline. 

(iii)  Wake-emitted  internal  waves. 

(iv)  Proper  orthogonal  decomposition  of  the  wake-radiated  internal  wave  field. 

(v)  Eulerian  mean  flows  during  the  reflection  of  an  internal  wave  beam  off  a  free-slip  surface. 

(vi)  Lagrangian  mean  flows  and  dispersion  in  a  reflecting  internal  wave  beam. 

(vii)  Subsurface  signature  of  the  wake-radiated  internal  wave  field. 


(i)  Mean  Flows  in  Large-Amplitude  Internal  Wave  Packets . 
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Figure  1:  Left  panel:  Schematic  of  problem  setup  where  a  vertically  localized  IGW  is  generated  in 
the  forcing  region  and  allowed  to  freely  propagate  downward.  In  the  critical  level  simulation 
(described  in  greater  detail  in  the  Ph.D.  thesis  of  A.M.  Abdilghanie)  ,  a  Gaussian  background  shear 
flow  is  positioned  near  the  bottom  of  the  computational  domain.  Right  panel  :  Two-dimensional 
contours  of  the  horizontal  velocity  field  after  five  wave  periods  in  a  fully  nonlinear  simulation  of  a 
horizontally  periodic  vertically  localized  wave  packet.  The  virtual  wave-generator  is  located  at 
z/2^2.5  where  a  significant  mean  flow  is  observed. 

2-D  internal  wave  beams  (IWB)  or  packets  are  often  used  as  idealized  models  for  internal  waves  radiated 
by  submerged  topography  in  the  ocean  and  the  internal  wave  field  emitted  by  a  localized  stratified 
turbulent  event  such  a  turbulent  wake.  Both  IWBs  and  internal  wave  packets,  nevertheless,  carry  some  of 
the  salient  features  of  the  above  geophysical  and  naval  internal  wave  fields,  including  their  ability  to 
support  significant  nonlinearity  and  dispersionA 
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Figure  2:  Comparison  of  the  mean  horizontal  flow  in  a  fully  nonlinear  simulation  and  linear 
simulations,  along  with  simulations  in  which  either  the  zero  Fourier  mode  is  removed  or  the 
proposed  mean  flow  absorbing  term  is  used,  for  a  horizontally-periodic  vertically-localized  wave 
packet. 
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Figure  3:  Vertical  profile  of  the  instantaneous  horizontal  velocity,  sampled  through  the  middle  of  a 
wave  packet,  compared  at  one  and  five  wave  periods  in  a  fully  nonlinear  simulation  of  a 
horizontally  periodic  vertically  localized  wave  packet  in  which  the  mean  flow  absorbing  term  is 
used. 


Our  very  first  effort  within  this  ONR-funded  project  aimed  to  study  the  interaction  of  a  finite-amplitude 
2-D,  horizontally  non-localized  internal  wave  packet,  with  a  background  shear  layer,  as  shown  in  the  left 
panel  of  Figure  1 .  Nevertheless,  we  found  that  our  wave  packet  underwent  significant  distortion  and 
decay  as  soon  as  it  was  generated  ;  it  fully  degenerated  far  from  the  source  (right  panel  of  the  same 
figure).  We  established  an  understanding  of  the  mechanism  underlying  this  mean  flow  formation  which 
was  found  to  lie  in  the  gradient  of  the  vertical  Reynolds  stress  field  of  the  wave  packet.  A  scaling  of  the 
strength  of  the  mean  flow  as  a  function  of  wave  amplitude,  stratification  characteristics  and  wave  packet 
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localization  was  provided.  This  scaling  was  then  used  to  construct  various  mean  flow  mitigation 
techniques  to  enable  a  robust  wave-amplitude  internal  wave  packet  propagate  sufficiently  far  from  its 
source  free  of  any  mean  flow  (figures  2  and  3).  This  wave  packet  ultimately  interacts  with  variations  in 
background  stratification  in  current.  Such  interactions  are  characterized  by  significant  nonlinearity  and 
can  provide  important  insight  into  the  remote  fate  of  wake-radiated  internal  waves.  More  detail  on  this 
study  may  be  found  in  the  article  by  Abdilghanie  and  Diamessis  ( Theor .  and  Comp.  Fluid  Dyn.  2012). 


(ii)  Internal  wave  beam-pycnocline  interaction 

The  driving  question  behind  this  particular  investigation  was  whether  a  significant  fraction  of  wake- 
radiated  internal  wave  energy  can  propagate  through  variations  in  the  ambient  stratification  profile.  This 
work  was  conducted  in  close  collaboration  with  Dr.  Scott  Wunsch  at  the  Johns  Hopkins  U.  Applied 
Physics  Laboratory.  Since  the  simulation  of  the  interaction  with  an  overlying  pycnocline  of  the  internal 
wave  field  radiated  by  an  actual/full  stratified  full  wake  would  be  a  prohibitively  costly  computational 
endeavor,  we  opted  to  focus  on  the  interaction  of  an  upward  propagating  internal  wave  beam  (IWB)  with 
a  pycnocline  above  which  also  sits  a  mixed  layer  (see  Figure  4  for  a  schematic  of  the  computational 
domain  and  an  example  stratification  profile. 
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Figure  4:  Left  panel:  Schematic  of  the  computational  domain  indicating  the  dominant  phenomena 
along  the  IWB  trajectory:  primary  and  secondary  reflection  and  refraction.  The  color  shading 
denotes  the  background  density  profile  extending  from  light  to  heavywith  increasing  depth.  Right 
panel:  Example  Brunt-Vaisala  frequency  profile  showing  the  two  fundamental  parameters  of  the 
problemunder  consideration:  Pycnocline  thickness,  /i,  and  the  ratio  r  =  Nmax/N\  of  maximum 
buoyancy  frequency  inside  the  pycnoclineto  its  value  in  the  lower  layer. 

Dr.  Wunsch’s  equivalent  experimental  studies  had  revealed  the  formation  of  high  amplitude  harmonics 
upon  encounter  of  the  IWB  with  the  pycnocline  ;  these  harmonic  waves  were  typically  trapped  inside  the 
pycnocline  and  tended  to  propagate  downstream  within  it.  To  this  end,  our  numerical  work  also  focused 
on  harmonic  formation.  24  simulations  were  conducted  examining  a  range  of  pycnocline  to  lower  layer 
buoyancy  frequency  ratios,  r  =  Nmdx/Nh  and  normalized  pycnocline  thicknesses  h/Ax,  where  Xx  is  the 
IWB’s  horizontal  wavelength.  Our  work  was  restricted  to  IWB  propagation  angles  of  45°  with  respect  to 
the  horizontal,  typical  of  what  is  observed  in  wake-radiated  internal  wave  fields  (see  also  section  (iii)  ), 
and  a  beam  with  relatively  high  steepness  of  2%  .  Such  an  IWB  has  sufficiently  high  amplitude  to 
undergo  a  variety  of  nonlinear  effects  within  the  pycnocline  but  is  not  strong  enough  to  be  subject  to 
breaking. 
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Figure  5:  Snapshots  of  instantaneous  JN2/N2(z)  at  steady  state  for  0.1.  Here  JN2  is  defined  as 
the  perturbation  to  the  buoyancy  frequency  associated  with  the  wave-induced  density 
perturbation.  Top,  middle  and  bottom  panel:  r  =  2,  r  =  6,  and  r  =  10,  respectively.  The  thin 
horizontal  solid  lines  delineate  the  interval  to  z$<  z<  Zo  +  h! 2,  where  zois  the  pycnocline  center. 
The  upper  part  of  each  panel  is  blank  because  of  the  near-zero  value  of  background  stratification 
requency  N2(z)  at  these  depths 
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Figure  6:  The  normalized  amplitude  of  the  primary  frequency  (left)  and  harmonic  mode  (right)  for 
simulations  with  h/Xx=  0.1  and  r  =  2  (top)  and  r  =  6  (bottom).  The  r  =  6  case  is  near  the  interfacial 
wave  resonance.  The  color  scale  is  logarithmic.  Data  are  sampled  at  a  coarser  spatial  resolution 
than  that  of  the  actual  simulations. 


Our  findings  are  best  described  in  the  context  of  the  two  limits  of  pycnocline  thickness  we  considered, 
h/Xx  =0.1  and  1,  thin  and  thick  pycnoclines  respectively.  A  visualization  of  the  flow  structure  for  h/Xx  =0.1 
is  shown  in  Figure  5.  Increasing  r  shows  a  reduced  capacity  for  the  incident  IWB  to  reflect  back  into  the 
lower  layer.  Distinct  modal  structures  are  observed  within  the  pycnocline  for  all  cases,  with  the  r=6  case 
being  the  most  efficient  in  enabling  the  propagation  of  these  structures,  linked  to  harmonic  frequencies, 
downstream  within  the  oycnocline.  At  r= 2,  the  pycnocline  buoyancy  frequency  can  only  support 
evanescent  harmonic  modes  and  at  r=10,  viscous  decay  of  the  strongly  refracting  beam  prevents  strong 
generation  of  harmonics  and  the  subsequent  propagation  within  the  pycnocline. 

The  thin  pycnocline  case  of  r= 6  was  found  to  be  optimal  in  terms  of  harmonic  generation  and 
propagation  as  further  supported  by  a  weakly  nonlinear  analysis  for  plane  internal  waves,  developed  by 
Dr.  Wunsch.  The  analysis  is  accounts  only  for  refraction  (and  not  reflection)  of  the  incident  IWB. 
Optimal  harmonic  generation  and  propagation  occur  when  the  wavenumber  and  frequency  of  the 
harmonic  wave  are  in  resonance  the  dispersion  relation  of  the  interfacial  mode  associated  with  the  thin 
pycnocline.  Evidence  in  this  direction  is  provided  by  Figure  6,  where  the  r= 6  shows  a  much  larger 
harmonic  energy  propagating  far  downstream  within  the  pycnocline.  Additional  evidence  of  this 
particular  resonance  is  provided  by  the  agreement  between  theory  and  direct  numerical  simulation  in 
terms  of  harmonic  amplitude  at  different  downstream  locations  within  the  pycnocline  (Figure  7). 
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Figure  7:  The  dimensionless  harmonic  amplitude  as  a  function  of  the  parameter  y  (a  measure 
of  pycnocline  strength  ;  y= 2  corresponds  to  r= 6)  for  the  numerical  simulations  with  different 
pycnocline  characteristics,  along  with  the  analytic  result  for  a  thin  pycnocline.  The  analytical 
result  from  our  weakly  nonlinear  theory  shown  as  a  dashed  black  line  (also  shown  are  results  from 
a  previous  theoretical  study  by  Thorpe  based  on  a  completely  different  mathematical  formulation). 
Filled  black  squares  show  the  numerical  result  for  a  pycnocline  thickness  h!i.x=  0.1  at  x/  Ax  =  6.5, 
while  filled  blue  circles  show  the  numerical  result  further  downstream,  at  x!  Xx  =  8. 
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Figure  8:  The  amplitude  of  JNz/NZifor  the  primary  (left)  and  harmonic  (right)  modes  for 
simulations  with  h/Ax=  1  and  r  =  6.  The  color  scale  is  logarithmic.  The  dotted  black  line  is  the  ray- 
trace  path  of  the  incident  wave;  the  solid  and  dotted  white  lines  show  ray-trace  paths  for  the 
harmonic  frequency  2co.  Visualization  window  is  focused  on  the  vicinity  of  the  pycnocline. 


In  the  case  of  thick  pycnoclines,  harmonic  generation  is  also  dictated  by  refraction  but  not  by  the  above 
resonance  argument.  As  shown  in  Figure  8,  the  strongly  refracting  IWB  undergoes  significant  viscous 
decay  upon  its  entry  into  the  pycnocline.  However,  non-negligible  harmonic  generation  is  observed 
during  the  initial  refraction.  Harmonic  waves  spawn  off  along  their  own  propagation  paths  and  a 
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complex  harmonic  structure  is  established  through  nonlinear  wave  interactions  further  downstream  into 
the  pycnocline.  Inviscid  simulations  are  needed  to  further  clarify  the  observed  flow  structure  and 
dynamics. 

More  detail  on  our  findings  in  this  topic  of  research  may  be  found  in  the  paper  by  Diamessis  et  al.  (. Dyn . 
Aim.  Ocean.  2014). 


(iii)  Wake-radiated  internal  waves 

Previous  ONR-funded  work  by  our  group,  enabled  by  our  in-house  designed  spectral  multidomain 
penalty  method  code,  revealed  that  at  sufficiently  high  body-based  Reynolds  numbers,  the  simulated 
stratified  wake  of  a  towed  sphere  followed  a  life-cycle  differently  than  what  is  commonly  perceived 
(Diamessis  et  al.  J.  Fluid  Mech.,  2011,  671:  pg.  52-95):  Buoyancy  does  not  suppress  turbulence  within 
one  buoyancy  period.  Instead,  the  horizontal  shear  layers  is  generates  are  sufficiently  fine  and  strong 
such  that  secondary  Kelvin-Helmholtz  instabilities  arise,  followed  by  energetic  secondary  turbulence, 
which  can  last  up  to  M^O(IOO). 

We  investigated  the  impact  of  this  prolonged  secondary  turbulence  on  internal  wave  radiation  by  the 
turbulent  wake.  Two  sphere-based  Reynolds  number  values,  /te=5><10  and  10*  were  considered  along 
with  a  range  of  internal  Froude  numbers,  Fr=  4,  16  and  64  (Computational  limitations  at  the  time 
prevented  us  from  reaching  Fr=  64  at  7?e=105).  Our  primary  findings  were: 

a)  At  high  Re ,  /he  persistence  time  of  internal  wave  radiation  is  also  prolonged,  with  significant 
high-frequency  waves  being  emitted  from  the  wake  for  as  times  as  late  as  Nt=$0  (Figure  9) 

b)  There  is  a  distinct  variation  in  the  phase-line-tilt  of  wake-radiated  waves  with  increasing  Re 
(Figure  10).  At  low  Re ,  typical  of  the  laboratory,  all  Fr  emit  waves  at  approximately  30°  from 
the  vertical,  an  angle  which,  according  to  theory,  suggests  that  wave-radiation  is  a  viscously 
driven  process.  It  is  at  this  angle  that  viscous  losses  are  minimized.  Low  Re  wakes  are  highly 
inefficient  in  depleting  the  mean  wake  flow  of  its  energy  through  wave  radiation.  In  contrast,  at 
high  Re ,  waves  are  radiated  at  angles  in  the  range  [45° 55°].  Such  angles  optimize  the  extraction 
of  momentum  from  the  wave-generating  source. 

c)  Consistent  with  the  above  observation,  we  found  that  high  Re  wakes  can  extract  up  to  a  total  of  6 
to  8%  of  the  initial  mean  wake  momentum  (Figure  11),  unlike  their  low  Re  counterparts.  Given 
that  naval  or  geophysical  wakes  operate  at  Re  4  to  5  orders  of  magnitude  higher  than  the  largest 
Re  considered  in  our  study,  it  is  likely  that  wave  emission  can  be  a  significant  source  of  energy 
that  may  be  comparable  to  viscous  dissipation  and  diapycnal  mixing. 

Additional  results,  not  shown  here  have  addressed  the  power  law  decay  rates  of  the  horizontal 
wavelength  of  the  most  energetic  wave  packets.  Furthermore,  we  have  identified  a  tendency  for  an 
increase  in  wave  steepness  with  increasing  Re  and  Fr  (see  also  section  (vii)  and  the  publication  by 
Abdilghanie  and  Diamessis,  J.  Fluid  Mech  2013). 
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Figure  9:  Isosurfaces  of  horizontal  divergence  (a  measure  of  the  internal  waves  radiated  by  the 
localized  turbulence)  for  a  stratified  turbulent  wake  of  a  towed-sphere  at  sphere-based  Reynolds 
number,  Re=  105,  and  Froude  number,  Fr=  4.  Shown  are  two  different  times,  Nt= 20  (left)  and  80 
(right).  Note  that  a  wake  at  the  same  value  of  Fr  but  at  a  /?e=5*103  typical  of  the  laboratory  has  no 
internal  waves  associated  with  it  at  the  later  time.  The  persistent  internal  wave  radiation  at  high 
Re  is  directly  connected  to  the  secondary  turbulence  inside  the  wake  core. 
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Figure  10:  Phase-line  tilt  angle,  0,  with  respect  to  the  vertical  of  wake-radiated  internal  waves  for 
simulated  stratified  turbulent  wakes  of  towed  spheres  at  two  different  body-based  Reynolds 
number  values,  i?e=5xl03  and  105(empty  and  full  symbols,  respectively).  Data  are  shown  for  three 
internal  Froude  number  values,  Fr= 4,  16  and  64.  The  low  Re  data  concentrate  around  0=30°, 
whereas  the  high  Re  data  initially  have  values  in  the  range  [45°  55°]. 
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Figure  11:  Left:  instantaneous  vertical  flux  of  internal  wave- transported  horizontal  momentum 
computed  at  vertically  offset  Ay-planes,  sufficiently  far  from  the  wake  edge,  for  the  cases  shown  in 
the  previous  figure.  Empty  and  full  symbols  correspond  to  low  and  high  Re ,  respectively  (circles 
and  triangles  are  Fr=  4  and  16  cases).  U0  is  the  mean  wake  centerline  velocity  at  Nt=2.  Right:  Time- 
integrated  momentum  transported  by  the  stratified  turbulent  wake  for  all  of  the  above  cases. 


(iv)  Proper  orthogonal  decomposition  of  wake-radiated  internal  waves 

A  subset  of  our  analysis  of  the  wake-radiated  internal  wave  field  explored  the  potential  of  Proper 
Orthogonal  Decomposition  (POD)  to  probe  into  the  generation  mechanism  of  internal  waves  by  a 
stratified  turbulent  wake  of  a  towed  sphere.  As  part  of  a  preliminary  effort,  in  collaboration  with  Profs. 
Roi  Gurka  of  Ben  Gurion  U.,  Israel  (now  at  Coastal  Carolina  U.  in  the  U.S.)  and  Alex  Liberzon  at  Tel 
Aviv  U.,  Israel,  we  isolated  our  7?e=5><10  and  Fr=  4  datasets  and  applied  POD  to  the  vorticity  and 
horizontal  divergence  components  on  a  series  of  snapshots  on  select  2-D  slices  on  xz  andyz  planes  within 
the  interval  10  <Nt<  35.  Our  findings  are  published  in  an  article  in  Physics  of  Fluids  (Diamessis  et  al. 
2010)  and  are  summarized  here: 

a)  Reconstruction  of  the  streamwise  vorticity,  cox ,  fields  on  select  yz  planes  from  individual  modal 
subranges  indicated  a  preference  of  each  subrange  for  a  particular  angle  of  wave  radiation  from 
the  horizontal  (figures  12  and  13).  In  particular,  lower  POD  modes  captured  layered-like  features 
inside  the  wake  core  and  more  horizontally-inclined  internal  wave  beams,  linked  presumably  to 
lower  frequencies  of  radiation.  Higher  POD  modes  showed  more  fine-scale  structure  inside  the 
wake  core,  presumably  originating  from  turbulence,  whereas  the  internal  waves  came  out  an 
angles  that  were  more  vertical. 

b)  Timeseries  along  radial  lines,  inclined  at  increasing  distance  from  the  wake  centerline,  provided 
preliminary  quantitative  support  for  the  above  qualitatively  observed  correlation  (Figure  14). 

More  recent  work  with  our  Israeli  colleagues,  funded  by  the  US-Israel  Bi-National  Science  foundation 
has  provided  further  probing  into  the  capacities  of  POD  to  separate  turbulence  and  internal  waves  in 
turbulence. 


10 


Figure  12:  Full  (panels  (a),  (c)  and  (e)  )  and  reconstructed  (panels  (b).  (d)  and  (f)  )  (ox  field  on  a 
yz  plane  located  at  the  middle  of  the  wake-centerline.  Reconstruction  is  performed  for  POD  modal 
subrange  1-3.  Top  panel:  M=15.  Middle  panel:  Nt= 25.  Bottom  panel:  Nt= 35.  The  dashed  black 
lines  are  inclined  at  45°  to  the  horizontal. 


11 


Figure  13:  Reconstructed  a)x  field  on  the  sam eyz  plane  as  in  the  previous  figure  located  at  the 
middle  of  the  wake-centerline.  Reconstruction  is  performed  for  POD  modal  subrange  4—6  (left 
column)  and  7-10  (right  column).  Top  panel:  A^IS.  Middle  panel:  Nt= 25.  Bottom  panel:  Nt=35, 
The  dashed  black  lines  are  inclined  at  45°  to  the  horizontal. 
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(a) 


(b) 


Figure  14  (previous  page):  Timeseries  of  the  full  and  reconstructed  cox  field  on  the  same yz  plane 
as  the  previous  figure  obtained  from  virtual  sensors  positioned  at  different  azimuthal  angles  on 
the  1st  quadrant.  Sub-figures  (a),  (b)  and  (c)  correspond  to  modal  subranges  1-3,  4-6  and  7-10, 
respectively.  In  each  sub-figure  the  top-left  panel  shows  time  series  at  wake  centerline.  The 
remaining  three  panels  of  each  subfigure  show  signals  at  specified  azimuthal  angles.  Blue  and  red 
lines  correspond  to  data  at  r=3D  and  4D,  respectively.  Black  lines  are  used  only  for  the  data  at  the 
wake  centerline  in  the  top  panel.  Solid  lines  represent  full  flow  field  and  nonsolid  lines  represent 
the  reconstructed  flow  field. 

(v)  Eulerian  mean  flows  in  a  subsurface  reflecting  internal  wave  beam 

As  a  first  step  in  the  second  leg  of  our  ONR-funded  research  and  as  part  of  an  ongoing  effort  to 
understand  the  fate  of  turbulence-radiated  internal  waves  far  from  their  source,  this  study  examined 
nonlinear  effects  in  the  reflection  of  a  2-D  finite-amplitude  internal  wave  beam  (IWB)  off  a  free-slip 
surface,  a  model  for  the  ocean  surface  on  a  calm  day  (Figure  15).  A  linear  stratification  was  considered, 
along  more  complex  ones.  We  focused  on  mean  flow  formation  but  also  investigated  the  formation  of 
frequency-harmonics  of  the  primary  beam  and  the  potential  for  parametric  subharmonic  instability  .  In  all 
these  efforts,  a  uniform  linear  stratification  was  considered.  A  subset  of  our  analysis  considered  similar 
effects  in  a  configuration  where  a  mixed  layer  laid  above  the  linear  stratification.  The  base  of  the  mixed 
layer  operated  in  an  equivalent  manner  as  the  free-slip  surface. 


Figure  15:  Schematic  of  a  2-D  IWB  reflecting  off  a  free-slip  surface. 
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Figure  16:  (a)  Streamwise  position  vs.  time  diagram  of  u-velocity  at  the  surface  off  which  the  IWB 

reflects,  (b)  Time-averaged  IWB-induced  mean  flow. 

The  primary  findings  of  this  study  were  the  following: 

a)  We  have  identified  the  structure  of  the  IWB-induced  Eulerian  mean  flow  over  the  corresponding 
region  of  the  water  column.  A  persistent  current,  between  10  and  50%,  of  the  horizontal  phase  speed 
of  the  beam  is  established  at  the  surface  (figures  16  (a)  and  (b) ). 

b)  The  amplitude  of  the  mean  Eulerian  surface  current  is  found  to  scale  linearly  with  the  IWB  amplitude 
A  and  tan where  6  is  the  phase-line  tilt  of  the  beam  w/r  to  the  vertical  (Figure  17).  Strongly 
hydrostatic  beams  (6  very  close  to  90°)  produce  the  strongest  mean  flows.  Very  good  agreement  with 
the  theoretical  predictions  of  Tabaei  et  al.  (J.  Fluid  Mech.  2005,  526:  pg  217-243)  is  found  for  all 
simulated  waves. 

c)  For  sufficiently  high-amplitude  beams,  a  parametric  subharmonic  instability  is  observed  after  a  long 
enough  time  (Figure  18).  Significant  energy  of  the  primary  IWB  is  spawned  off  to  secondary  beams 
of  lower  frequency,  which  are  more  horizontally  inclined.  Frequencies  and  wavenumbers  of  the 
primary  and  secondary  beams  are  found  to  match  theoretically  prescribed  resonance  conditions  very 
closely. 
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Figure  17:  Left  panel:  Magnitude  of  IWB-induced  mean  current  at  the  free  surface,  normalized  by 
the  magnitude  of  the  maximum  horizontal  velocity  inside  the  primary  IWB,  as  a  function  of  beam 
amplitudes!  and  angle  to  the  vertical  0 .  Numerical  computed  mean  currents  vs.  the  theoretical 
estimates  of  Tabaei  et  al.  (2005)  as  a  function  of  angle  0 . 
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Figure  18:  (a)  Streamwise  position  vs.  time  diagram  of  instantaneous  horizontal  velocity  at  the 
surface  for  a  very  high  amplitude  beam  (st=6.8%).  Time  and  position  are  normalized  with  IWB 
period  T  and  horizontal  wavelength  Xx,  respectively.  Parametric  subharmonic  instability  (PSI) 
manifests  itself  in  the  form  of  more  horizontal  like  isophase  lines  at  t/7>50.  (b)  x:z-contour  plot  of 
horizontal  velocity  at  late  times  for  a  PSI-prone  IWB.  The  subharmonic  waves  are  more 
horizontally  inclined  and  are  easily  spotted  spawning  off  to  the  left  of  the  primary  beam. 
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(vi)  Lagrangian  mean  flows  in  a  surface  reflecting  internal  wave  beam 


Stokes 


Figure  19:  Leading  order  Lagrangian  mean  drift  (top)  vs.  Eulerian  mean  flow  for  the  reflecting 
IWB  in  figure  16.  Both  sets  of  contour  values  are  normalized  by  U02Icx ,  where  U0  is  the  maximum 
horizontal  velocity  in  the  incident  wave  and  cx  is  its  horizontal  phase  speed.  The  two  different  types 
of  mean  flow  cancel  each  other  out  to  almost  machine  precision. 

As  we  continue  to  probe  into  the  subsurface  signature  of  the  wake-radiated  internal  wave  field,  we  have 
further  extended  our  study  of  the  reflection  of  an  individual  IWB  at  a  free-slip  surface  to  explore  for  the 
Lagrangian  mean  flows  in  the  reflection  zone.  Understanding  such  Lagrangian  effects  is  imperative  to 
quantify  the  true  mean  flow  (as  its  Eulerian  component  does  not  suffice  in  this  regard).  Furthermore  it 
provides  insight  into  the  internal  wave-driven  Lagrangian  drift  of  surface-particulates. 

Borrowing  elements  of  surface  wave  theory,  we  have  found  that  the  Stokes  drift  fully  cancels  out  its 
Eulerian  counterpart  in  the  reflection  region  (figure  19).  Any  resulting  Lagrangian  particulate  drift  should 
be  of  second  (higher)  order  than  each  of  these  mean  flow  components.  Application  of  a  high-accuracy 
(Fourier-Legendre  in  space,  4th  order  Runge-Kutta  in  time)  particle  tracking  technique  to  the  flow  fields 
used  in  Zhou  and  Diamessis  (2013),  confirms  this  claim  (figure  20).  However,  the  dispersion  becomes 
significant  when  PSI  ensues  (figure  21). 
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Figure  20:  Tracer  particle  orbits  for  five  wave  periods  within  a  reflecting  IWB.  No  noticeable 
mean  drifts  are  observed.  Any  higher-order  Lagrangian  drift  is  very  weak  for  a  stable  reflecting 
IWB  at  this  particular  wave  set-up. 


Figure  21:  Tracer  particle  displacements  at  the  initial  time  (top)  and  after  one  wave  period 
(bottom)  for  the  PSI-prone  reflecting  IWB  of  figure  16.  Any  higher-order  Lagrangian  drift  is  very 
weak  for  a  stable  reflecting  IWB. 


18 


(vii)  Subsurface  signature  of  the  wake-radiated  internal  wave  field . 


This  is  ongoing  work  which  we  expect  to  wrap  up  in  late  summer/early  fall  2014  and  should  give  rise  to 
one  more  manuscript.  We  consider  the  same  cases  as  in  section  (iii),  though  in  deeper  and  wider  domains 
and  augmented  with  a  Fr= 64  simulation  at  the  higher  Reynolds  number.  Focusing  on  a  uniform 
stratification,  which  extends  up  to  the  surface,  we  aim  to  identify  the  time  (in  buoyancy  units)  of  when 
the  strongest  strains  emerge  at  the  surface  above  the  submerged  wake.  Furthermore,  we  seek  to 
parameterize  these  strains  as  a  function  of  Re  and  Fr.  In  addition,  we  are  exploring  the  variation  in 
observed  surface  wavelengths  and  their  orientation  across  all  runs  while  we  are  also  investigating  surface 
signatures  of  interest  to  remote  sensing  experts.  The  latter  effort  is  in  close  collaboration  with  Prof.  Bill 
Philpot  at  CEE-Comell. 

The  first  step  in  analyzing  the  (sub)surface  signature  of  the  wake-radiated  internal  wave  field  has 
consisted  of  calculating  the  surface  strain  field,  as  represented  by  the  horizontal  divergence  Az . 
Examination  of  the  7?e=5x10^  dataset  indicates  that  the  peak  strain  rate,  across  all  Fr ,  scales  as  Fr]  2 
whereas  the  corresponding  arrival  time  t/(D/U)  scales  as  Fr2  3  (figure  22)  The  equivalent  analysis  at 
7?e=105  is  still  underway.  A  first  finding  is  that  the  internal  wave  induced  strains  persist  for  a  longer  time 
and  reach  a  peak  value  that  is  0(10)  larger  than  their  low  Re  counterpart. 

A  key  component  of  ongoing  work  aims  to  quantify  the  characteristic  time  and  length  scales  of 
the  near-surface  flow  field  through  1-D  and  2-D  continuous  wavelet  transforms.  On  account  of 
the  significant  vertical  distance  separating  the  edge  of  the  wake  and  the  free  surface,  significant 
dispersion  (variation  in  frequency)  is  observed  with  lateral  offset  at  the  surface  (figure  23). 
Finally,  figure  24  shows  a  sample  of  2-D  wavelet  analysis  of  the  surface  strain  field  with  a 
distinct  resonance  occurring  at  a  particular  horizontal  wavescale. 

What  is  most  intriguing  is  the  structure  of  the  wake-radiated  internal  wave  field  at  7te=10s  and  Fr= 64 
(figure  25).  Along  with  the  above  dispersion,  the  wake  energy  released  into  the  internal  wave  field  is 
sufficiently  large  that  the  emitted  waves  are  so  steep  that  they  are  very  close  to  breaking.  The 
implications  of  such  potential  breaking  for  surface  signatures  will  be  explored  in  the  near  future. 

The  surface  signature  of  wake-radiated  internal  waves  may  be  remotely  detected  through  various  means, 
one  of  which  is  the  characterization  of  slicks  of  surface  surfactants.  By  solving  for  a  pure  2-D  advection 
equation  on  the  surface  we  can  observe  the  evolution  of  a  passive  scalar  (with  an  initially  uniform 
distribution)  as  affected  by  the  wave-induced  surface  velocity  field.  The  slick  pattern  mimics  that  of  the 
reflecting  internal  waves.  Most  importantly  though,  the  enrichment  ratio  (instantaneous  over  initial 
concentration)  of  the  passive  scalar  field  is  found  to  exceed  the  critical  value  of  1 .05  at  Re=105  (figure 
26)  which  implies  enhanced  remote  detectability  of  the  surface  signature  of  the  waves  due  to  strong 
visible  contrast 
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Figure  22:  Surface-averaged  internal  wave-induced  strain  at  the  free  surface  for  the  /?e=5xl03 
datasets  at  three  different  Fr  values  (left).  Right  panel  indicates  a  preference  for  a  preferred  Fr- 
scaling  for  strain  magnitude  and  time  of  observation  of  maximum  strain  at  the  surface. 
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Figure  23:  Top:  Streamwise  position  vs.  time  (in  buoyancy  units)  diagrams  of  horizontal 
divergence  at  the  free  surface  for  the  /?e=5xl03and  Fr= 4  simulation  at  three  different  lateral 
offsets.  Bottom:  Corresponding  jc-averaged  period.  Notice  the  steady  change  in  wave  period  with 
lateral  offset,  an  indicator  of  frequency  dispersion  in  the  original  wake-radiated  internal  wave 
field. 
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Figure  24:  Top:  Snapshot  of  surface  horizontal  divergence  for  the  same  simulation  as  in  the 
previous  figure.  Only  the  positive  j-half  of  the  free  surface  is  shown.  Bottom:  2-D  transect 
from  the  continuous  wavelet  transform  modulus  cube  at  they/D  offset  corresponding  to  a 
white  line  in  the  top  panel.  Red  regions  of  strong  resonance  correspond  to  individual  wave 
packets  of  horizontal  wavelength  Xh/D. 


Figure  25:  Isopycnal  contours  from  an  yz  slice  from  the  /?e=105and  Fr= 64  simulation  at 
M=45.  Notice  the  very  steep,  nearly  overturned,  isopycnals  in  the  radiated  internal  wave 
field  at  non-negligible  distances  from  the  wake  core.  Whether  these  waves  break  is  a  topic 
of  current  investigation. 
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Figure  26:  Top  panel:  Horizontal  divergence  field  at  two  different  times  for  the  ^?e=10sand 
Fr= 4  simulation.  Bottom  panel:  Corresponding  concentration  fields  resulting  from  a  2-D 
advection  simulation  with  an  initial  condition  of  uniform  scalar  concentration  (values 
normalized  with  initial  concentration,  a  quantity  defined  as  enrichment  ration  in  remote 
sensing).  Notice  how  the  enrichment  ration  exceeds  the  critical  value  of  1.05  at  late  times, 
which  suggests  enhanced  likelihood  of  remote  detectability. 
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Ongoing/Future  Work 


Current  work  is  focused  on  analyzing  the  full  subsurface  signature  of  wake-radiated  waves  for  the  six 
pairs  of  points  in  the  ( Re,Fr )  parameter  space  identified  in  section  (vii)  of  this  proposal.  Completion  of 
this  effort  is  anticipated  by  late  summer  2014,  to  be  followed  by  the  write-up  of  a  relevant  publication. 
Through  the  assistance  of  DoD-PETTT  we  are  working  on  profiling  and  optimizing  the  performance  of 
our  parallel  spectral  multidomain  penalty  method  solver  for  the  incompressible  Navier-Stokes  equations. 
Enabled  by  a  DoD  Frontier  project,  spearheaded  by  Prof.  Steve  de  Bruyn  Kops  at  U.  Mass.  Amherst,  we 
are  engaged  in  regular  discussions  on  stratified  turbulence  and  have  begun  implementing  direct 
numerical  simulations  (with  Prof,  de  Bruyn  Kops)  initialized  with  our  own  implicit  L.E.S.  data.  In 
parallel,  through  Frontier  computational  resources,  we  have  initiated  three  massively  parallel  stratified 
wake  simulations  at  /?e=4><105  and  Fr=  4,  16  and  64.  Run  completion  is  expected  in  late  October  2014 
with  the  associated  analysis  completed  by  late  winter  2015.  The  analysis  will  aim  to  quantify  the  de¬ 
scaling  of  properties  of  the  mean  wake  and  the  radiated  internal  wave  field. 
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